Thermal diode: Rectification of heat flux 
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By coupling two nonlinear one dimensional lattices, we demonstrate a thermal diode model that 
works in a wide range of system parameters. We provide numerical and analytical evidence for 
the underlying mechanism which allows heat flux in one direction while the system acts like an 
insulator when the temperature gradient is reversed. The possible experimental realization in nano 
scale systems is briefly discussed. 
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The understanding of the underlying dynamical mech- 
anisms which determines the macroscopic laws of heat 
conduction is a long standing task of non-equilibrium sta- 
tistical mechanics Recent years have witnessed some 
important progress in this direction even though a satis- 
factory understanding is, so far, unavailable. For exam- 
ple, it has been surmized that for a one dimensional lat- 
tice, momentum conservation leads, in general, to anoma- 
lous heat conduction Q, with heat conductivity k which 
diverges with the system size L as L 2//5 0]. If the trans- 
verse motion is considered, then k ~ L 1 ^ 3 4]. A connec- 
tion between the anomalous conductivity and anomalous 
diffusion has been also established^. Moreover, after 
two decades of debates, it is now clear that exponential 
local instability is not a necessary condition for the va- 
lidity of Fourier law[(|. 

A better understanding of the mechanism of heat con- 
duction may also lead to potentially interesting applica- 
tions based on the possibility to control the heat flow. 
Indeed, recently, a model for a thermal rectifier has been 
proposedQin which the rectifying effect is obtained by 
acting on the parameters which control the nonlinearity 
of the lattice. Although this model is far away from a re- 
alistic implementation, nevertheless it opens the possibil- 
ity to propose thermal devices which may have practical 
relevance. 

In this Letter we demonstrate, via computer simula- 
tions, the possibility to build a thermal diode by cou- 
pling two nonlinear lattices. Indeed our model allows 
the heat flow from one end to the other, but it inhibits 
the flow in the opposite direction. The model works in 
a wide range of system parameters and the underlying 
microscopic mechanism is different from the one in ref|jj 
even though it is based on the same general idea. In the 
present case the ratio of the heat currents in the two op- 
posite directions is 100 times larger than that in modelpj 
and it is therefore sufficiently large to encourage experi- 
mental work. 

Our system consists of two segments of nonlinear lat- 



tices coupled together by a harmonic spring with con- 
Each segment is described by the 



stant strength fcj„ t 
Hamiltonian: 
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In fact, Eq. (QJ is the Hamiltonian of the Frenkel- 
Kontorova (FK) model which is known to have normal 
heat conduction 0. 

For simplicity we set the mass of the particles and the 
lattice constant m = a = 1. Thus the adjustable param- 
eters are (k L ,k int ,k R ,V L ,V R ,T L ,T R ), where the letter 
L/R indicates the left/right segment and T LR is the tem- 
perature of the left/right heat bath. In order to reduce 
the number of adjustable parameters, we set V R — XVl, 
k R = Xk L , T L = T (l + A),T R = T (l - A) and, unless 
otherwise stated, we fix Vl — 5, fcz, = 1 so that the ad- 
justable parameters are reduced to four, (A, A, ki n t, To). 
Notice that when A > 0, the left bath is at higher tem- 
perature and vice versa when A < 0. 

In our numerical simulations we use fixed boundary 
conditions and the N particle chain is coupled, at the 
two ends, with heat baths at temperatures Tj, and T R 
respectively. We use Langevin heat baths and we have 
checked that our results do not depend on the particu- 
lar heat bath realization (e.g. Nose-Hoover heat baths). 
We then integrate the differential equations of motion by 
using the 5th-order Runge-Kutta algorithm as described 
in We compute the temperature profile inside the 
system, i.e. the local temperature at site n defined as 
T n = m(x n 2 ), where ( ) stands for temporal average, 
and the local heat flux J n = k(x' n (x n — x„_i)) [lj. The 
simulations are performed long enough to allow the sys- 
tem to reach a non-equilibrium steady state where the 
local heat flux is constant along the chain. 

In Fig. ^ we plot the heat current J versus A for 
different temperatures To. It is clearly seen that when 
A > the heat current increases with A, while in the 
region A < the heat current is almost zero, i.e. the 
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FIG. 1: Heat current J versus the dimensionless temperature 
difference A for different values of To- Here the total number 
of particles N = 100, k int = 0.05, A = 0.2. The lines are 
drawn to guide the eye. 



system behaves as a thermal insulator. The results in 
Fig. ^shc-w that our model has the rectifying effect in 
a wide range of temperatures. The rectifying efficiency 
depends on temperature as well as on other parameters 
as described below: 

ki n t effect. The interface elastic constant ki n t is a very 
important parameter as it plays the role of coupling left 
and right lattices. By adjusting this parameter one can 
control the heat flow. Indeed, once the parameters of the 
two lattices are fixed, then the smaller the coupling the 
smaller is the heat current through the system. 

In order to describe quantitatively the rectifier effi- 
ciency we introduce the ratio, \J+/J—\ where J+ is the 
heat current (from left to right) when the bath at higher 
temperature T + is at the left end of the chain and </_ is 
the heat current (from right to left) when the left end of 
the chain is in contact with the bath at lower temperature 
T_. 

FigEK a) shows that, by varying ki n t> the rectifier effi- 
ciency | J+/J_| changes from about two times at k^t = 
0.5 to more than 100 times for ki nt < 0.01. More im- 
portantly, this figure shows that the rectifying effect is 
very significant in a wide range of ki nt . In the inset of 
this figure we show | J± | versus ki„t for a system of eight 
particles. It is seen that, for small ki n t, numerical data 
follow the dependence J± oc kf nt over almost two orders 
of magnitude. Therefore the ratio | J+| / |J_| is a con- 
stant independent of kint- 

Fig. Otb) shows the temperature profile for different 
ki n t ■ There exists a large temperature jump at the inter- 
face and this jump is much larger when the high temper- 
ature bath is at the right end. In this case there is very 
small temperature gradient inside each lattice and, as a 
consequence, the heat current is almost vanishing. 

Effect of the lattice parameter A. As it is known 
in the FK model the elastic constant and the strength of 
the on-site potential can be scaled to a single parameter. 
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FIG. 2: (a) Heat current |J±| versus the interface elastic con- 
stant kint for N = 100, A = 0.2, T+ = 0.105, T_ = 0.035. 
The solid circles are the current J+ while open circles are 
the current J_. Inset in (a) is \J±\ versus kint for a system 
of eight particles. The solid lines in the figure and in the 
inset have slope equal two. (b) The temperature profile for 
kint = 0.01, 0.05, 0.2 for the same parameters of (a) 



Therefore it is sufficient to study the properties of the 
system (1) as a function of the single parameter A. 

In Fig. |21 we show the current J± versus A for two 
different interface constants h m t — 0.05, and kint = 0.2. 
This figure clearly show that in a wide range of parame- 
ters our model has a quite good rectifying efficiency, i.e. 
\J+/J- \ ~ 100. 

Rectifying mechanism. To understand the underlying 
rectifying mechanism, let's start from the energy spec- 
trum of the interface particles. Fig. 0] shows the phonon 
spectra of the left and right interface particles at dif- 
ferent temperature when the two lattices are decoupled 
(hnt = 0). 

The match/mismatch of the energy spectra of the two 
interface particles controls the heat current. It is clearly 
seen from Fig|2|that, if the left end is in contact with the 
high temperature bath Tl , and the right end with the low 
temperature bath Tr (< Tl), then the phonon spectra of 
the two particles at interface overlap in a large range of 
frequencies, thus the heat current can easily go through 
the system from the left end to the right end. However, 
if the left end is at lower temperature Tl and the right 
end is at higher temperature Tr (> Tl), then the phonon 
spectrum of the right interface particle is mainly in the 
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FIG. 3: Heat current versus A for kint = 0.2 (circles), and 
kint = 0.05 (triangles). The solid symbols refer to Tl = 
T+,T R = T-, while empty symbols refer to T L = T-,Tr = 
T+. Here T+ = 0.105, T_ = 0.035 and N = 100. 



low frequency part, while the left interface particle is in 
the high frequency part. Then there is almost no overlap 
in phonon frequency, and the heat current can hardly 
go through from right to left, and the system behaves 
as an insulator. Why the left and right particles at the 
interface have so different phonon spectra? This can be 
understood from the following analysis in different tem- 
perature regimes. 

(1) Low temperature limit. At low temperature, the 
particle is confined in the valley of the on-site potential. 
By linearizing the equation of motion one can easily ob- 
tain the frequency band 7]: 



V < lu < VV + Ak. 



(2) 



For the case of FigH with T=0.01 (left) and T=0.002 
(right), this corresponds to 0.36 < ui/2ir < 0.48 for the 
left particle and to 0.16 < u)/2-k < 0.21 for the right 
particle. 

As the temperature is increased, the interparticle po- 
tential kx 2 /2 becomes more and more important until 
a critical value T cr w V/(2tt) 2 is reached (we take the 
Boltzman constant equal unity), when the kinetic energy 
is large enough to overcome the on-site potential barrier. 
At this point low frequency appears and this happens at 
the critical temperatures T cr = 0.13 for V = 5 (left), and 
T cr = 0.025 for 7 = 1 (right). This is in quite good 
agreement with the data of Fig. 0] 

(2) High temperature limit. In the high temperature 
limit the on-site potential can be neglected, the system 
is close to two coupled harmonic chains, and the phonon 
band is [lflj . 



< lj < 2Vk, 



(3) 



which gives < lj/2tt < 0.32 for the left particle and 
< lo/2h < 0.14 for the right particle, again in good 
agreement with Fig. 0] 
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FIG. 4: Spectra of the two particles at the interface for dif- 
ferent temperatures at ki n t = 0. (a) particle at the left side 
of the interface, (b) particle at the right side of the interface. 
Here A = 0.2, N = 100. 



In fact, in order to optimize the rectifying effect, one 
should avoid the overlapping of the phonon bands in the 
low temperature limit (Eq|5J and that in the high tem- 
perature limit (Eq|3J for each segment of the system. Ac- 
cording to the above estimates, one should have V > 4fc, 
which is satisfied for the case of Fig. 0] 

The analysis presented here allows also to understand 
why in Fig^there is a region of negative conductivity for 
T = 0.09 and -0.5 < A < -0.2. Indeed in this interval 
it is seen that by increasing the temperature difference 
between the two heat baths, the current decreases! As 
explained above, the reason for such unusual behavior 
is due to the fact that by decreasing the temperature of 
the left bath, the phonon spectra of the two interface 
particles separate one from the other. 

It is worth pointing out that it is the temperature de- 
pendence of the phonon band that makes the rectifying 
effect possible. This dependence is due to the nonlinear- 
ity of the potential and therefore it should be possible to 
observe the rectifying effect in any nonlinear lattice. The 
difference between different choices of the on-site poten- 
tial is quantitative rather than qualitative. For example, 
we have observed the rectifying effect in the tfi 4 model 
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FIG. 5: The finite size effect of the rectifying efficiency. Here 
we plot |J+| (solid circles) and |J_| (open circles) versus N. 
Parameters are A = 0.2, k int = 0.05, T+ = 0.105, and T_ = 
0.035. 



wherein the on-site potential is kx 2 /2 + /3a; 4 /4, but the 
rectifying efficiency is not as significant as it is in the FK 
model. The main reason for the very large rectifying effi- 
ciency in the FK model is the finite height of the on-site 
potential. 

Finite size effect. The results shown above are for a 
system with total number of particles N = 100. Since the 
rectifying mechanism in our model is due to the coupling 
between two dissimilar lattices it is reasonable to expect 
that the number of nonlinear oscillators will definitely 
affect the rectifying efficiency. This is shown in Figure 
where we plot | J + | and | J_ | versus the system size N 
with T+ = 0.105, and T- = 0.035. 

The decrease of the rectifying efficiency by increas- 
ing N (Fi^SJ and k int (Fi^2J can be understood in the 
following way. If we denote by R the sum of thermal 
resistances of the two segments, then the total resis- 
tance of the system is R± = R + R int , where Rf nt is 
interface resistance when A > and A < 0, respec- 
tively. As J± = (Tl—Tr)/R±, the rectifying efficiency is 
\J+/J-\ w (R+Rf nt )/(R+Rt nt )- Notice that R increases 
with N while R int are N independent, they decrease with 



increasing kmt, and R^ nt is much smaller than R int . 



It is 



therefore clear that the efficiency decreases with increas- 
ing ki n t (Fi^EJ) at fixed N and by increasing N at fixed 
k^t (Fi©. 

For the data of Fig|S]we have R~ nt 3> R, and therefore 
J_ is a small constant that is almost independent of N, 
while J + ~ 1/N. As N is increased over a critical value 
N c where \J+\ |J_|, the rectifying effect vanishes. For 
the case of Fig. [5] this happens at N c ~ 10 5 . However, 
we should point out that the value of N c depends on ki n t 
and by decreasing ki n t we may increase N c . 

Finally, we would like to discuss the possibility of an 
experimental realization of our diode model. In our 
calculations, we take dimensionless units. For a typi- 
cal atom, the dimensionless temperature Tq used in this 



paper is related to the real temperature T r by T r ~ 
10 2 - 10 3 T 8]. Therefore T - 10° is of the order of 
room temperature. If we consider the lattice distance of 
1 A, then a lattice of 100 particles is about 10 nanome- 
ters long, which is of a size scale that can be fabricated 
in nowaday's laboratory. Another problem for real ex- 
periments is that the on-site potential is due to the in- 
teraction with a substrate. In our numerical simulations, 
we neglect the heat current in the substrate. This may 
affect the rectifying efficiency of the diode. Therefore, 
one needs to choose a substrate with substantially low 
thermal conductivity. 

In summary, we have devised a thermal diode by us- 
ing two coupled nonlinear lattices. Our model exhibits a 
very significant rectifying effect in a very wide range of 
system parameters. The underlying mechanism has been 
discussed as well as a possible realization of the model in 
a nano scale system. 
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